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This thesis is part of 
Postgraduate School to 
elctroraacnetic noise on tie 
station. The project's ooj 
signals in the 0.05 to 10 3 
the importance of obtain! 
nci se on the sea floor and 



the long term effort of the Naval 
collect and interpret EL? 
sea floor and from a nearby land 
ective is to study ana interpret 
z frequency ranee. It emphasizes 
ng measurements of getffiaar.se Lc 
on land over a oeriod of several 



years . 

The specific obectives of this thesis are twofold; to 
calibrate, install, and operate tie land data collection 
system; and to develops a computer program to analyze the 
data recorded to provide easily utilized interpretive 
^ ^ 0 



A. TIME VARIATIONS 0? THE GE03AGN: 



The earth's magnet 
for navigation at sea 
These common lses are 



c field has 
and for surv 
de pendent o 



field. The static field can be gener 
magnetic dipole at a 11.5 degree inc 
rotation of the earth. it shows onl 
tions in direction. Its n agni t ude / 
earths surface but usually is in t 
70, 000 nano-Teslas (nT) . At Monterey 
field has a value of 48,600 nT. 



oeen used for centuries 
ey direction on land . 
a the s t a z i c magnetic 
ally characterized as a 
lination zo zhe axis of 
y minor diurnal varia- 
aries greazly over the 
he range of 40,000 to 
, California the static 



It has only been 
tier, that the underiy 
become well known. 



sou rces . 



The sources 



iace the discovery of magnetic indue - 
ng variations In the static field have 
These variations are due to many 
include, but are not limited "C, solar 
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activities, lunar daily variations 



, magnetic storms, thun- 



derstorms in the earth's atmos 
detonation of nuclear devices, 
to this is the interactioi of t 
sun with the earth's magnet osph 
in duration, amplitude, an 1 the 
Within the frequency range of 
variations are characterized as 
Gecmacnetic micr cpulsa t ions 
tions in the earth's magnetic 
0.2 seconds to 1 0 minutes. They 



phere, ear 


t hquak es , and the 


One o: the 


major contributors 


he charged 


particles from the 


ere. These 


events vary widely 


f requ ency 


range they effect. 



this system the most common 
micr o o ulsatior.s. 

are iefi.ced to be fluctua- 
ted! which have periods from 
vary in amplitude from 0.001 



iiT tc as high as a few tens or nl. They are attributed to 
the interaction of groups of charged particles propagating 
in the earth's magnetosphere as hyirn m a gr.etic waves. These 
waves inturn create the scarce currents, (primarily in the 
icnsphere), for the micropulsation field. Micropulsar ions 



are divided into two subc 
nuous . The abr-eviatior s 

con ri rue us and pi for pu 
difference between the Own 
have a semi- sin usoi dal r.at 
show irregularities in boo 
period ranges and average a 
listed in Table I, [ Hef. 1 ] 
The overwhelming major 
frequency range are due i 
"Schumann" resonances, will 
of 6-8 Hz. These resonances 
in the concentric spneri 
surface and the lower regio 




is that the dc micropulsations 
ure while tie pi micropulsations 
h amplitude and frequency. The 
mpiitude for th»se pulsations are 



i ty of the pulsations in this 
o solar activity. En addition, 
also be apparent ia the vicinity 
are due to lightening transients 
cal cavity between the earth's 
n of the ionosphere, *Rsf. 2]. 
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r ABLE I 



Geomagnetic 1i cropulsncl on Classes 



Notation Period(5eo) Amplitude (nT) 



pci 


0. 2-5 
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5-10 


0.1-1 


pc3 


10-45 


0.1-1 


DC4 


45-150 


0.1-1 


DC 5 


150- 600 


1-10 


bi 1 


1-40 


0.01-0. 1 


?i 2 


40-150 


1-5 



SYSTEII DESCRIPTION 




Figure 1.1 Sensing Coil Dimensions 



ID 



The sensor system currently employed consists of a thrs* 
coil land system deployed at La Mesa Village housing area, 
approximately 2 km from tie Naval Postgraduate School data 
collection center, and a two coil sea system deployed on the 
bottom of Monterey 3ay. Its distance to the Naval 
Postgraduate School (NPS> varied from 2 to 4 km. The 



original intention of the project v as to deploy 



:he 



.an: 



system approimately 22 miles from the Naval Postgraduate 
Scnocl at Chew’s Ridge. This proved impossible due to data 



See A p o e n d i x 



^ ~ * ior.a 1 



transmission difficulties, 
ir.f crma~icn on this site. 

Each sensing coil is a coa tin ao a sly wound coil antenna 
made up cf approximately 5460 turns of 13 gauge copper wire. 
The average sensing area of each coil is 0.0825 square 
meters and its dimensions are as depicted in Figure 1.1. 
The sea system has its coils mounted orthogonally, usually 
in the horizontal plane. The land system coils are affixed 
to a wooden mount that holds them ii an orthogonal configu- 
ration; see Figure 1.2. The directional orientation of the 
sea system cannot be determined when deployed. The land 
system is oriented wita the X- coil in the magnetic ncrta- 
sou th direction , Y-coil in the magnetic east- west direction 
and the Z-coil in the vertical. 
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Figure 1.2 Sensor Wanting Blocic 

The signal coni i tin ners , PCM ooard = vi transmitter a:s all 
housed in a small wooden 3 a il ding and powered by D . 3 . power 
supplies. This building is the closesr source of alternating 
current . 

After the data has been encoded oy the PCM, it is tran- 
smitted by VHP radio data link to tie Naval Postgraduate 
School for recording. 

C. DATA COLLECTION AND PROCESSING 

Land and sea data are simultaneously received and 
recorded on separate channels of aa analog tape recorder. 
Each tape holds 60 to 120 minutes of data. These tapes must 
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Figure 1.3 Land Data Collection System 



be decoded and re-recorded cr.no dig 
tha main- frame computer an the Chum 
digital tape contains 93 ninutss of 
of these digital napes is then 
ccnputer program, (See Cnapter 3, 
the plcns of the various interpretive 
3 for rh e derails of the retorting 
cording on diaitai tape. 



inal tape suitable for 
Computer Center. Each 
Lind or S9i d 2 . t 2 . • Einh 
onocrsssd through th? 
ration CO/ to provide 
factors. See Reference 
, decoding and re-re- 
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II. SYSTEM CALIBRATION 

A. INTRODUCTION 

The triaxial coil magn etometer system is diagramed in 
Figure 2.1. The lata sensed by this magnetometer system is 
given by V x (t) , Vy ( t) and l z (ti ; the output voltages gener- 
ated by the input magnetio fields. Tie system transfer func- 
tion must be determined so that the magnetic field 
com pone n.ts can be removed from the output voltage. In the 
time domain the relationsaip between the output voltage and 
the magnetic field is a convolution integral, 

v; (*-) s f h; (t-t'l o; (flit- (1) 

where i=x,y,z. It is more convenient to express Sgn. (1) in 
the frequency domain. 

V; {«) =H; (ui)3 / M . (2) 

The transfer function Hi can easily be determined in 

terms of the ratio of output tc an input test signal cr by 
analysis of the measurement system’s properties. 

The next section discusses the components of tne system 
and attempts to develop a system model. The last section 
describes the calibration measurements and the experimental 
determination of the transfer function. 

3. SYSTEM DESCRIPTION 

The individual transfer function for the sensing coil 
is: 



V c (c«) = j**»K3 (*u) 



(3) 
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Figure 2.1 Systea Coiponent Tor-figuration 



where V c (u») 


is the spectram of tie 


the coil and 


K is a gain factor. 


The r.ex* 


co m po n er.f in the s y s t e m 


pass filter. 


The preamplifier intro 



voltage directly out of 

is the preamplifier low 
duces a gait of 1000 so 



the voltage signal can be sent to remote instrumentation. 
The lew pass filter is a five pole Chebyshev with a 3i3 
cutoff at 20 Hz. The filter is primarily used to remove 50 
Hz signals from the system and to prefiiter the data before 
if is digitized. An additional pole is introduced because of 
a capacitive coupling between the second and third stages of 
the filter. Figure 2.2 shows a pole-zero diagram of the 
preamplifier. The transfer function takes the form: 



H (j^ =K/ ( jw+a) (jw + b) (jo/fc-jd) ( j<*/ + r j 3 1 (j^+e-jf! (j«/+a + jf) ('4| 
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where K is the gain cf tie ore ampli f i er and the coeficisnts 
a - f are determined from the pole-zero diagram. 

The preconditioner nas additional gain cf approximately 
1030 and has zen ar diodas to limit tie input voltage to the 
PCS board. Typically the cutoff is around 7 to 3 volts. 

The magnitude and phase response of the preamplifier- 
filter and pr eco nditicr.a r are shown in Figures 2.3 through 
2.3. Th c input test signal generating these curves had flat 
magnitude and phase character as snowr. in Figures 2.9. The 
figures then represent the actual transfer functions of the 
preamplifier-filter preconditioner system for the x-y-z 



channels of land and ocean sensor systems. The figures show 
that the chan ne Is are well matched in magnitude and phase. 
The phase response is linear and virtually identical in 
every case. This is an important point since the phase par 1 - 
cf the transfer function is not removed from the output 
voltage signal. However, since the data analysis is in terms 
of power spectral densities ar.d coherences, this fact has no 
bearinc cn the results. 



The total transfer function of 
of Zgn. (3) and (U) . 3ased on these 
the resulting transfer function is 
linear at the low frequency end. 



e system is the product 
ormulas and the figures 
predicted to be fairly 



C. DETERMINATION 0? TH3 SYSTEM TRANSFER FUNCTION 



The magnitude of the system trans 
mined experimentally through tie use 
source cf sinusoidal current and a s 
experimental set up is illustrated 
turn Helmholtz coil with a 0. SI me 
meter separation between coils was u 
bra ted magnetic field to tie sensing 



fer function was deter- 
of a Helmholtz coil, a 
pact rum analyzer. The 
in Figure 2.10. A one 
ter radius and a 0.51 
sed tc provide a cali- 
co il. The sensing coil 
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Pigure 2.2 Pole-Zero Diagra* for 3od = l 13-10A Aaplifi?r 
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Figure 2.3 Response of the Land X System 
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Figure 2.4 Response of the Lani Y System 
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Figure 2.5 Response of the Lard Z Systa 
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Figure 2.6 Response of the Ocean X System 



21 



X. IS. 823 Y.ra.7M Y. 25. 287 

TRANS IU 3 I 




Magnitude R -sp;n s : - 



X, 22. 557 Yi 27. 7 72 Y, Ci. #45 

"RANS #A* 3 




Figure 2.7 Response of the Dcsan Y System 
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Figure 2.8 Response of -he Dcean Z System 
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Figure 2*9 Response of the Input Test Signal 
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was placed concentric to the Helmholtz coil, and centered 
vertically, to ensure a uniform magnetic field about the 
sensing coil. The magnetic field in this region of the 
Helmholtz coil is directed along its axis and is given by 

[Ref. 4 1 : 

B N 1/5 (5) 

to terms on the crier of (L/R), where L is the linear dimen- 
sion of the region, !1 is the number of turns, R is the 
Helmholtz ceil radius and I is me driving current of the 
coil. The driving current was produced by a Wavetek signal 
generator, Model Number 142. The current was measured by 
reading the voltage across a 996 ohm resistor in series with 
the Helmholtz coil. A spectrum analyzer (HP-3592A) was used 



to measure 



a -u , 



system outpu J 



to avoid 



any SDurious signals 
generated by harmonics. 3y repeating the measurements at a 
number cf different frequencies and drive currents, the 
system transfer function could be measured. 

The actual calibration was conducted ac the La Mesa 
village site using “he setup illustrated in Figure 2.10. 
The Helmholtz ceil, sensing coil and preamplifier were 
located approximately 40 meters from the shack which housed 
the remainder of the experiment. S' it: h the Helmholtz coil 

disconnected the background magnetic noise of this location 
was repeatedly examined to ensure accurate readings. Figure 
2.11 shows a typical background noise plot from the speerrum 
analyzer (the transfer function has not been removed). The 
distinct features are Schumann resonances, otherwise the 
background was low enough and stable enough to adequately 
calibrate the system. 3nly at the Schumann frequencies was 
the background subtracted from the measured output voltage. 
Calibration runs were conducted from 0.05 Hz to 20.0 Hz with 
applied calibrated magnetic field of 0.02, 0.2 and 1.0 
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Figure 2.10 System Calibration Experiment 



nancteslas. Figure 
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calibration curves 


of the land 
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ocean magnetometer 


systems. The true 


transfer fa net i 


o?. is 


represented by the 


0.02 nT field. The 


larger applied 


fields 


show saturation of 



the system at the protective zener diodes in the precor.di- 
tior.er. At the low frequency end all the systems are closely 
matched. At 1 Hz an output voltage of 1 v indicates a 
measured magnetic field of 0.4 nT. A 1 nT field is saturated 
above 2 Hz. 

The computation of the equations, representing the 



magnitude 


of 


the transfer 


functions, requires 


fitting the 


data to 


the 


theoretical tr 


ansfer function of 


the previous 


ssc t , i on . 


Th 


is would be de 


sirable if the phase 


character of 
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the sys-em transfer function nasi be estimated. However, 
since the primary area of interest is in the lower frequency 
region, and since the phase response is not needed if the 
various channels maybe assamed identical, it was decided to 
model the transfer functions in terms of linear segments. 
The equations derived in this manner have a maximum error of 
1^ in the range of 0 to 1 0 Hz and 5^ in the 10 to 20 Hz 
range for the maonitude of the transfer funtion. 

In the actual computation each curve was broken into 
convenient segments that could be represented by a straight 
line. The range of the s egm en.ts were 0 to 5 Hz, 5 to 7.5 Hz, 
7.5 to 10 Hz, 10 to 15 \Z and 15 to 20 Hz. The straight 
line parameters were complied for =*ach segment to yield a 
series of five equations to model each sensing ooil and 
amplification sys-em. See Table El for the individual 
transfer functions. 



27 




Figure 2.11 Background Oaomagnetic Noise f (Hz) 05/20/82 
15 13 hrs. 
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Figure 2.15 Ocean System X-Ooil Calibration 
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Figure 2. 16 Ocsaa System Y-Soil Calibration 
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r ABLE II 



Transfer Function Compu-.-r ?rc<jr-.n 



CC ? L = 1 ,N 
F F C- t-R £G (L ) 

I F ( F R 0 • IE.25 . )GC TC 1 
xMii-n(ii/ 2 a. 

GC T r 3 

1 IF (F3Q. Lc .15. ) GC TC 2 

XX (L >=XX (l ) / ( 1C5.5-3.I4“FF*J) 

YY = YY(L)/ (131. 32-7. 533* P*C) 

Z ( U -l ( L J / ( I 7 7. 2c- 7. A3 4r F FQ) 

GC TC 3 

2 IF (FRO. L E . 10 . ) GC TC 3 

XX{U = XX(L»/(5.«5e*FRQ-3C.97> 

Y Y ( L > = Y Y ( L ) / {7 . 1 1 fc + P RQ-3 9 . 09 ) 
Z(U = Z (U/<o. 4 S*FPC- 32. 35) 

C-C TO 3 

3 IF (F?q. l E .7. 5 > GC TC A 

>><U = X*(L>/(3.AC2*FRQ-6. 2 1 > 
YY{l>=YY(l)/(4.2;~*FPC-:0.35> 

Z ( L ) =Z ( L ) / (A .OAA^FCC-7 .3 9 ) 

GC TC 3 

4 ie 3.LE.5. )GC TC 5 
XXa)=XML)/(2.fc?ll*PPJ+0.14 O fc7) 
YY { L ) = Y Y { t_ ) / ( 3 .G1Z^FRQ-1. c 5) 

Z { l )=/. ( L ) / (3 . 1 8A» FFC-1 ) 

GC TC 9 

5 IF(t-Rj.Lt.i. )GC 'L 6 

>X(U*XX<L)/(2.*>2 1I«FRO0.L4Sf.7) 

7 YY(U*YY(L)/( 2.7C2*F'<gi 
Z ( L I = Z ( L )/ {2.C2*FRC > 

GC TC a 

6 XX (L ) = X x (L )/ < 2 . 72 ♦FRO 

GC TO 7 

S CCNTPJir 
<7 CCNTIM'JE 
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DATA EVALUATION 

A. POWER SPECTRAL DENSITI AND ALIASING 

The first step in the analysis of the data collected was 
the determination of the power spectral densities (PSD’s) 
for each channel of the land ar.d sea lata. This was accom- 
plished through the use of the computer program developed by 
J. H. Schweiger [Ref. 5]. After minor modification of this 
program, approximately twenty hours of real time data was 
evaluated. Samples of the PSD's are presented in Figures 3.1 

S.0 — . .. » 

The overall value cf the PSD’s is in agreement with 
previous work conducted in this area. The PSD's have a 
decline cf approximately - 6 d3 per octave from 0 to 5 Hz 
where they tend to decrease their rate of decline. Night 
PSD’s tended to be 3 to 7 iB below tae daytime level. This 
coincides with the findings cf HcDevi v t and Hcmar. (1980), L-l. 
W. Beard (1981) and A. C. Eraser-Smith and J. L. Buxton 
(1975), [Ref. 6]. 

One feature of the PSD's that differed from previous 
PSD's is the spike at 4 HZ in the land data. This was found 
to be "aliasing" of the 53 Hz power grid component by the 
PC.d sampler. This experiment was conducted with a sample 
rate of 32 Hz; which, by che sampling theorem, should give a 
clear view of the frequency range of 0 to 10 Hz, when 
coupled with the 20 Hz cut-off filter in the ampli fier , [ Ref . 
7]. The 60 Hz component cf the geomagnetic field was 
measured by F. W. Slayton ( 1979) and found to be 4 (r.T) /Hz, 
which is 10 higher than the surrounding data [Ref. 8]. The 
20 Hz filter, built into the amplification system, would 
normally eliminate frequencies in this range but because of 
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Figure 3.1 P3D, X-Coil, Land Data 
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Figure 3.2 PSD, T-Coil, Land Dana 
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Pigara 3.3 PSD, Z-Coil, Land Dara 
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Figure 3.4 ?SD r X-Cail, Sea Da"a 
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Figure 3.5 PSD, T-CoiL, Sea Da-a 
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the great differences in magnitude this does not happen. 
When the 32 to 64 Hz range is convolved onto the 0 to 32 Hz 
range the 60 Hz component will appear at 4 Hz. This destroys 
the usefulness of the data in a small region about this 
point . 

Another point of interest is illustrated by the vertical 
component of the sea data recorded on September 16, 1982, 

Figure 3.5. During this experiment the Y-coil was oriented 



in a vertical rather than 


th e 


more nor 


mal 


hori 


zontdl dir e c ~ 


tior. . Examination of the 




shows 13 


to 


nS 

to 

r— 


B decrease in 


th 9 PSD. This is in a arse 


n enr 


with : h e 


theory 


than, because 



of the ocean -a ir interface, geomagnetic waves propagate 
vertically in the sea, thus eliminating this component. 



3. STOKES PARAMETERS AND THEIR INTERPRETATION 



After tne cc 
Seotetnber of 193 2, 
a prime concern, 
computer program 
Spectral Densities 



lection of the 
the problem of how 
During the latter 
was developed to 
(PSD » s| by J. H. 3 



data during June to 
"to i^ reverse it b ^ c a in ~ 
earn of this r e r i o d a 
dstsoinins ths ? c w 9 r 
hweiaer, October 1982, 



[Her. 5]. 



This program provided the basic structure to read the 
data points off of magnetic data tape and, through the use 
of a Fast Fourier Transform (FFT) routine and transfer func- 
tions, to compute the geomagnetic fluctuations in the 
frequency domain. This program provided the basic format 
which was expanded and modified to produce other parameters. 
The derivations of the ecu a* ions for these parameters and 
their interpretation comprise the remainder of rhis section. 

The first parameter to be discussed is the coherence 
between the different channels on land and sea. Coherence 
gives a measure of the "singleness" cf the source at the 
various frequencies. "Singleness" in this context refers to 
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Figure 3.6 PSD, Vertical Component of See Data 



U3 



what percentage cf the intensity at a frequency can be said 
to come from a single source. To compute the coherence one 
first calculates the cross power spectral densities between 
the channels. For a given pair of channels X and Y, the 
cross power spectral density(S ) is defined: 

N 

V f J )=v/>X(fj ) Y (fj ) 
n= 1 

In this equation N is the number of blocks o 
through the FFT end averaged; X { f j ) and Y (f ;) a 
along the appropriate axis. The co herence (t^ r ) 
as: 

S*<Sj ) 

<*J)= • 

(S x (f j ) S 7 (f j ) ) 

For the land data coherences were also determined for tne 
Y-Z plane and the Z-X plane. Figures 3.7 and 3.3 show 

typical land and sea coherences, with all coils oriented in 
the horizontal plane. h'ote that on the land the coherence 
stays very close to its maximum value of one until approxi- 
mately 7-8 Hz and then starts to decrease while in the case 
cf the sea data there is more structure at the lower 
frequencies. This type of behavior is typical of most of the 
data, although the exact frequency of the decline varies. 

In addition ■‘-.o coherence, Stokes parameters may be 

computed from the auto-spectral densities and cross-spectral 
densities. They can be used to interpret the geomagnetic 
waves on a bi-Linear basis, a bi-Circuiar basis or an ellip- 
tical basis. 

Although other notations are utilized for the Stokes 
parameters, the most common notation; sO, si, s2, s3; will 
be employed here [Ref. 9]. The general form of *he Stokes 
parameters are: 
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Figure 3.7 Land Coherence, X Vs Y Channels 
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Fiqure 3.8 Sea Soherenca, X Ys Y Chanaal. 
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- 01.5 - 01.0 - 0(15 00.0 0(15 01.0 01.5 



sO=<a1 2 >+<a2 i > 
s1=<a1 2 >-<a2 2 > 
s2=2<a1 a2 zos S > 
s3=2<a1 a2 sin<£ > 



where al ana. a2 are the instantaneous amplitude of the 
orthogonal components of the magnetic vector and <T is their 
phase difference. In terms of the auto- and the cross-spec- 
tral densities the Stokes oarameters oar. be written: 



s0= 2 (5, +S y ) 



s2=4Re (S* y ) 
s3=4 Im ( Syy) 

It must be noted that the Stokes 
originally to describe monochrome 
applied to cur problem after the 



parameters were developed 
io light and can only be 
ata has beer, transformed 



to the frequency domain. 



The bi _ Linear basis irtertretatior of the 
tecs can be seen by examining a superposition 
on.ai linearly polarized waves [Ref. 10 1- 
Stokes parameters in this basis yields: 

s0= 1/2 (Ex<-Ey ) 
s1=1/2 (B*-bJ) 
s2= 3^ ByCOS (0 r -& r ) 
s3= B x By sin (9y ) 



Stokes para me - 
o f two or t hog- 
computing the 



Another basis useful for interpretation is the 
bi-Circular basis. This basis is represented by the superpo- 
sition of two counter-rotating circularly polarized waves. 
In this basis the Stokes parameters are: 



sO = B^. 
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Figure 3.9 Basis for Interpretation 

s1=23 + E_ cd s -j>_ ) 
s2 = 23^_ B_ sin ( -<p_ ) 
s3=3^- B± 

k final basis of interpretation is the elliptical basis. 
In this case the basis is modeled by th<= major and minor 
axis vectors and the angle formed between the X-axis and the 
major axis of the ellipse. From this basis the Stores param- 
eters are: 



s0=1/2 

s1 = 1 /2 (B^-5 M ) cos 2 ^ 
s2=1/2 (3 M -B„) sin 2 *c 
so = 3 W 
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The important point to none from these different basis 
secs is than for a given w a ve form (linear, elliptical, or 
cicular) the Stokes parameters, as computed from the spec- 
tral densities, are the same in any basis in which they are 
computed. This can be easily visualized by taking a circu- 
larly cr linearly polarized wave and determining the value 
of the Stokes parameters in each basis. It seems than these 
interpretations are most easily visualized in the elliptical 
basis. Consequently, in addition, to the Stokes parameters 
the ancle sxr was computed to make t he interpretation easier. 
This angle was calculated from the Stokes parameters = 1 and 
s2 by the equation: 

ok =1/2 arctan (s2 /s 1 ) 

Examples of the results of the calculations of the Stokes 
parameters and the angle <X for land ar.d sea data are 
displayed in Figures 3.10 to 3.19. Sots that si, s2, and s3 
are normalized by dividing by sO. In addition, for the land 
data, the same graphs were produced for the y - Z plane and 
the Z-X plane. 

The use. of Stokes parameters is based cn monochromatic 
and polarized waves. lo determine the applicability to our 
data of the Stokes parameters the degree of polarization, 
must be found. For a perfectly polarized monochromatic wave 

_a z .2. 

sO — s 1 +s2 + s3 . 



If the wave is partially polarized then sO becomes greater 
than the sum of the squares of the other three [RSF.9]. To 
determine how close to the polarized case the data lies, the 
degree of polarization (P) can be calculated. It is defined 
by. 



p= 



, ,2 - 2 . .2 
(si +s2 +s3 ) 



s3 
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Figure 3.10 Land SO, x-Y Plan- 
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Figare 3.11 Land S1/S3, X-Y Plans 
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Figure 3.12 Land S2/S0, X-7 Plane 
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Figure 3.13 Land S3/S0, X-Y Planr 
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Figure 3.14 Lind Angleoc, x-Y Plans 
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Pigure 3.15 Sea SO, X-Y Plane 
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Pigura 3.15 Saa SI/SO, X-Y Plana 
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Figure 3.17 S?a S2/S0, X-Y Plan* 
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Figure 3-18 Sea S3/S0, 1-1 Plane 
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Figura 3.19 Sea inglacC, X-Y Plaaa 
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If P is one then the wave is perfectly polarized and any 
analysis based on the Stokes parameters will be completely 
accurate. 

For values of P less than one the analysis is corres- 
pondingly less valid. However, factors that depend only on 
si, s2, and s3, such as the angle and the ellipticity, are 
net effected by low values of P. In such regions the noise 
level is much higher and the data is rot easily interpreted. 

The land a r.d sea lata vary in their values of ? 
<3 eg ending on th a rime of day and oner unknown factors. In 
general they both shew a high degree of polarization through 
3 Hz while above this frequency a marked drop off occurs. 
If was also noted that in the land data the degree of polar- 
ization is virtually identical in all three planes. This is 
aoparentlv due to -the fact that the polarization was not 
oriented along any of the three sensing axis and therefore 
the proportion of the total intensity that r s su It e d from 
polarized waves remained relatively constant. (See Fraures 
3.2 0-3. 23) 



Cne final 


too 


I cf 


int er prefa 


zi o n 
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ellipticity of 
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(2) . The 
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of icitv 


the major and 
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Specifically 2 


i s 


define 


d by: 










tan 


2=B m /B M 






This can also 


b e 


wr if ta 


r. in terms 


0 h 


the Sto 



as: 



2=1/2 arcsin (s3 / ( s 1 Z +s2 2 + s 3 2 ) x ) 



where E varies from 77V U to - . The interpretation of 2 is 

illuminated by considering its dependence on s3 and the 
elliptical basis. I f 2 is zero then the minor axis is zero 
which means the polarization of the wave is linear, at E= 
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Figure 3.21 Degree of Polarization, Land, T-Z Plane 
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Figure 3.22 Degree of Polarization, Land, z-X Plane 
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Figurs 3.23 



Degree of Pol* rizaf ion , 5aa 
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Figure 3.24 Elli priri ty, Land, X-Y Plane 
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Figure 3.2 5 
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the polarization is circalar and is rotating in a right- 
handed sense, and at E=- ^/4 the wavs is again circular, but 
with a left-handed rotation. In between these values of E, 
the polarization of the wave is an ellipse with the appro- 
priate handedness of the polarization depending on whether S 
is positive or negative. The smaller the absolute value of 
the ellipticity the narrower or more linear the ellipse 
becom es. 

From Figures 3.2 4 and 3.25, the form of these waves for 
both lard and sea data is generally a narrow ellipse with 
varying left and right handedness in the areas of high 
polarization while appearing to be more elliptical outside 
these areas [Ref. 9]. 

The final result of this investigation is that through 
the computation of the Stokes parameters the nature of the 
geomagnetic wave can be determined. Specifically, tc find 
the type of polarization o~ the wave, one computes the 
ellipticity and the orientation angle o< . The degree of 
polarization is also needed tc es t a biish the degree of 
validity of the interpretation. Naval Postgraduate School 
Technical Reoort 5NPS 61-33-005, contains additional plots 
of all ~he factors described in this section. 

C. COMPUTER PROGRAM DEV ELD PEMENT AND DATA PROBLEMS 

As noted earlier the ability to process large amounts of 
data in a timely manner was first put to use in this project 
in September 1982. The original program has been expanded 
to compute all of the factors listed in the previous 
section. It consists of a main program, a EFT subroutine, 
and a read subroutine. As noted before the FFI subroutine 
transfers the data into the frequency domain and the read 
subroutine reads the PCM data from the digital computer tape 
and places it in the input array. All changes and additions 
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made were in the main program, therefore only the main 
program is included in Appendix A. This program, Landxyz, is 
used to evaluate only the land data. The version of the 
program used to evaluate the sea data is identical to 
Landxyz with the following exceptions: sea system transfer 

functions are used; because the sea system has only two 
coils, cr.ly one plane of interpretive factors is computed. 

Although the basic program proved adequate to find the 



power spect 


ral densities, the growth 


i of the pregram 


* n 


include the 


ether calculations created 


several problems. 


The 


Jcb Control 


Language (JCL) orovides 


the guidance to 


the 



operating system of the com purer to tell in hew :o run batch 
jebs. Specifically what input and output devices to use and 
hew much memory to allocate • t :o the program during each step 
of its execution. (In the program the JCL can be identified 
by the //*s at the start of each line cf this code.) During 
the expansion of the program the amounts of memory in the 
LINK, GO and PLOT steps in the execution cf the program have 
all expanded. The current JCL is capable of handling up to 
two hours of data at 32 Hz or one hour of data at 64 Hz. 

Sea and land data are initially recorded on analog tape 
and then decoded and re-recorded onto digital tape in the 
Near Shore Laboratory of the Naval Postgraduate School. This 
procedure is discussed in [Ref. 3]. The output of this 
re-recording is approximately ninty minutes of real time 
data, from either the sea or land site, on a digital tape 
with a density of 800 bits per inch (BPI) . The storage 
density proved a great handicap in the rapid expansion of 
the program. The Church Computer Tenter has fifteen tape 
drives, but only one of these is capable of handling a 800 
BPI tape. The competition for this one tape drive led to 
processing times of up to twenty-seven hours. This delay 
made any programming error extremely costly in terms of time 
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lost. An attempt was male to transfer selected tapes into 
the memory of the computer, but it was unsuccessful due to 
minor recording errors on the digital tapes. These 
recording errors were overcome by the read subroutine in 
order to be able to work with the data. Time did not permit 
the presuit of this avenue of improvement but at best only a 
few tapes could be stored this way because of the vast 
amount of raw data. 

The computer procram was validated through the use of a 
test tape. The test taoe was recorded with a sinusoidal 
input at known frequency and voltage level as measured by a 
spectrum analyser. When this tape was processed by the 
Landxyz program, the correct frequency and amplitude were 
produced, [Ref. 11]. 

A noise floor tape was produced by replacing the sensing 
coils with 120 ohm resistors. This tape was then evaluated 
by the Landxyz program to produce the PSD of the noise floor 
of the system. Figure 3.25 shows this noise floor superim- 
posed cn a land PSD plot. 
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Figure 3.25 Noise Floor PSD 



67 



IV, 



CONCLUSIONS AND 3ECDMHENDATIONS 



A. CONCLUSIONS 

As stated earlier part of the objective of this thesis 
was tc calibrate, install and operate the land portion of 
the Naval Postgraduate School’s geomagnetic data collection 
system. This was generally successfully acomplished with the 
exception of degradation of data produced by using the mors 
noisy La Mesa site rather than the more remote Chew’s Ridge 
site. Exactly how much of an improvement in the quality of 
the data will result from the move to Chew’s Ridge is diffi- 
cult tc predict. It could possibly diminish the 50 Hz compo- 
nent so that data could be obtained in the region of u Hz. 
Overall some 70 hours of land data was collected during this 
research, of which approximately twenty hours has been eval- 
uated to date. 

For the thesis’s second objective, the production of a 
computer program and procedures to interpret the data, the 
expansion of the basic program to compute the coherences, 
Stokes parameters, degree of polarization, elliptical orien- 
taion angle, and the ellipticity would appear to have satis- 
fied the requirement. The major drawback to this system is 
the limited amount of data that can be processed in a 
reasonable time. At most the current programs are capable of 



evaluating on 


ly 5-6 hour 


s of land 


and sea data 


in a twenty- 


four period. 


This is 


pri marily 


iu? to •’■ape 


drive avaii- 


ability, but 


the reguir 


ement of t 


he man-hours 


f or decoding. 


re- recording 


on digital 


tape and 


filing these 


tapes at the 



computer center cannot be discounted. The data processing 
system has reached the stage where large amounts of data can 
be evaluated but it is not yet ready to operate on a conti- 
nuous basis. 
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RECOMMENDATIONS 



The first recommendation is that the data transmission 
problem be solved so that the Ihew’s ridge site can be oper- 
ated on a continuous basis. This has two advantages; in the 
improvement in the quality of the sampling environment and 
in giving a contiuous data stream that can be recorded at 
appropriate intervals to provide a long term data base from 
which low frequency trends could be observed. 

Improvement could also be made in "he area of data 
processing. Specifically, a new method of data handling is 
needed. Several alternatives exist, the first being the 
acquisition of the protests own decoder and digital tape 
drive tc provide a digital tape of a higher density than can 
now be made. This would make it possible for a two to eight 
fold reduction in the r.umoer of tapes necessary to hold the 
data and would also give faster turn — around time during 
processing due tc the increased numoer of tape drives avai- 
lable an the computer center. Other data handling schemes 
are also feasible if the problem o? how tc read data from 
the digital tapes into the mass storage system of the 
computer can be overcome. This would make if possible tc run 
programs in other than the ’’batch" mods and avoid the use of 
tape drives to gain access tc the data. 
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APPENDIX A 



COMP J TER PROGRAM 



//LANDXYZ JOB ( 1 57 3 , 0 1 29 ) , ' POGUE SMC 1 573 • f CL A SS = F 
//* MAIN OR G=NPG V Ml . 1 573P,LI NES=(99| 

//♦FORMAT PR r DDN AMS = PLOT . 3 Y SVECTR, DS3 T = LOCAL 
// EXEC FRTXCLG ?, REGION. G0= 3000K 
//FORT. SYSIN DD ♦ 

INTEGERS IN (16) 

C ARRAY 'IN' IS USED IN READING DATA FROM TAPE 

COH?LEX*8 XX ( 3 1 92) , YY (3192) ,2 (3 192) 

COMPLEXES SXY (8192) , SYZ (3 192) , SZX (8192) 

C THE COMPLEXES ARRAYS ARE USED TO ORDER INPUT 

C DATA AND INITIALLY REPRESENTS VOLT AGE- TIME 

C SERIES INFORMATION 

DIMENSION TIME (319 2) , FREQ (8192) , WORK (16 384) , FHQ 2 (31 92) 
DIMENSION ZX ( 3 1 92) ,ZY (8192) ,22(3192) 

DIMENSION SOXY (819 2) , SO Y2 (8 192) , SOZX (8 19 2) 

DIMENSION S1XY (819 2) , SI YZ (31921 , S1ZX (8192) 

DIMENSION S2XY (819 2) , S2 YZ (819 2) , S2ZX (819 2) 

DIMENSION S3XY (8192) , S3 YZ (8192) , S3ZX (8192) 

DIMENSION COXY (3 1 9 2) , CO YZ ( 8 1 9 2) ,COZX (8192) 

DIMENSION AXY (8192) , A YZ (3192) ,AZX (8192) 

DIMENSION PXY (8192) , ? YZ (8192) ,?ZX (3192) 

DIMENSION DEXY (819 2) , DEYZ (8192) ,D SZX (8192) 

C THE 'Z» ARRAYS REPRESENT FREQUENCY DOMAIN (FFT) 

C MAGNITUDE DATA AND ARE EVENTUALLY CONVERTED TO 

C POWER SPECTRAL DENSITY INFORMATION. COXY IS COHERENCE 

C MEASURED X TO Y. AX I IS THE ANGLE BETWEEN THE MAJOR 

C AXIS OF THE ELLIPSE AND THE X AXIS. PXY IS THE DEGREE 

C OF POLARIZATION IN THE X-Y PLANE. 

C DEXY IS THE SLLIPTICITY IN THE X-Y PLANS. 
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c 

c 



c 

c 

c 

c 

c 



INTEGER*'-* IT 3 ( 12) /12* 0/ 

REAL*4 RT3 (28) /29*0. 0/ 

SEAL ALAB ( 11) /» X’ , • Y' , ' Z' ,'C', 'SO', 'SI', '52', 'S3', 'A' 

,' P* , 'E'/ 

REAL* 8 TITLE (12) 

EQDIV AL ENCE (TITLE ( 1* , RT3(5) ) 

ARRAYS ' ITB' , 'RT3' , ' ALAB ', AND 'TITLE' ARE USED IN 
THE VERTS ATEC PLOTTER OUTPUT. 

DATA XX , YY/16 3 84* (0. , 0.)/ 

DATA Z/S1 9 2* (0 . ,0. ) / 

DATA ZX,ZY/16384*0./ 

DATA COXY, SOXY/16384* 0./ 

DATA 32XY, S3XY/163 34 * 0 . / 

DATA SI XY, S1YZ/1633U* 0./ 

DATA COYZ , SOYZ/16384 * 0./ 

DATA S2YZ, S3 Y 2/ 16 3 84 * 0 . / 

DATA COZX , SOZX/16384* 0. / 

DATA S2ZX , S3ZX/16384* 0. / 

DATA TZ,?REQ/16364*9. / 

DATA AXY, AYZ/1 6384*3. / 

DATA AZX, S 1ZX/ 16384*0 ./ 

DATA PXY, PYZ/1 6384*0. / 

DATA PZX,DEXT/16384*0 ./ 

DATA DEZX, DEYZ/16384* 0./ 

PI 2=6 .2831 853 

THE FOLLOWING SECTION READS THE DECODED PCM DATA OFF 
THE DIGITAL TAPE AND FOR EACH CHANNEL SKIPS EVERY 
OTHER SAMPLE POINT. 

I5EC=3 
ITL=ISEC*3 2 
ITL=300 

DO 55 JJ=1 ,ITL 

CALL RD <20 ,IN ,200,IREC,IRR) 

55 CONTINUE 



71 



IFRAME=81 9 2 
NR= 16 

FNR=FLOAT (NR) 

DO 70 L 1= 1 f NR 

C THE DO LOO? ENDING WITH STATEMENT 70 ENABLES THE PP.O- 

C GRAM TO PROCESS A LARGE AMOUNT OF DATA BY REPEATING 

C THE PROCESS IN BLOCKS. THE DATA POINTS FROM EACH RUN 

C THROUGH THE DO LOOP ARE ADDED TOGETHER AND EVENTUALLY 

C AVERAGED EY THE THE NUM3ER OF TIME THROUGH THE LOO?. 

C ’NR’ REPRESENTS THE NUMBER OF DATA SEQUENCES TO 5E 

C AVERAGED 

C 1 SEQUENCE CURRENTLY EQUALS 3192 DATA POINTS FOR EACH 

C CHANNEL OR 256 SECONDS AT A SAMPLE RATE 3? 32 CPS. 

DC 60 JJ= 1 ,IFR A ME 

CALL RD (20 ,IN, 1003 ,IR EC,IRR) 

XX (JJ) =IN ( 2) 

YY (JJ) =IN ( 3) 

Z (JJ) =IN(4) 

60 CONTINUE 

WRITE (6 ,2 0 0) IRR r IREC 

2 00 FORMAT ( 10 X , ’ IR P= ' , 15, 5X,' IREC = ' ,16,/) 

C THIS SECTION GENERATES THE TIME AND FREQUENCY 

C ARRAYS AND NORMALIZES THE INPUT PCM DATA TO VOLTAGE 

C IN PREPARATION FOR FAST FOURIER TRANSFORM TO THE 

C FREQUENCY DOMAIN. 

N = 8 1 9 2 
FN=FLOAT ( N ) 

DELT AT= 1 . / 32. 

T=FN*DELT A T 
DELTA?=1. 0/T 
DO 20 J=1 , N 

TIME (J) =D ELI AT * FLO AT ( J) 

FREQ (J) =D ELTAF *FLOAT ( J) 

XX (J) = (XX (J) -2 045. 5) * 10. /2 04 5. 5 
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XX (J) =REA L (XX (J) ) 

YY (J) = (YY ( J) -2 045. 5| * 10. /2 04 5. 5 
YY (J) =REAL (YY ( J) ) 

Z (J) = (Z (J) -204 5.5) *1 0 ./2045. 5 
Z ( J) = REAL (Z(J) ) 

C ’XX’ IS THE X-COIL DATA, 'YY 1 15 THE Y-COIL DATA, 

20 CONTINUE 

WRITS (6,9 0 0) X, Y,Z? 

201 CONTINUE 

DO 2 1 J = 1 , N 

FRQ2 ( J) =A L OG 1 0 (FREQ ( J ) ) 

21 CONTINUE 

C THE NEXT 3 STATEMENTS PERFORM AN FFT ON IHE INPUT 

C TIME SERIES DATA. 

CALL FOURT (XX, N, 1,-1, 0, WORK) 

CALL FOURT (TY,N,1 , -1 , 0, WORK) 

CALL FOURT (Z , N , 1,-1, 0 ,NORK) 

C THIS SLOCK OF STATEMENTS APPLE THE SYSTEM (VOLTAGE TO 

C E-FTELD) TRANSFER FUNCTION TO THE TRANSFORMED 

C FREQUENCY DOMAIN DATA. THIS 3L0CX ENDS AT 9. 

C THE TRANSFER FUNCTION CONVERTS VOLTS TO NANOTESLAS. 

C tjig ***** 

C THIS TRANSFER FUNCTION IS GOOD FOR MAGNITUDE ONLY. 

DO 9 L= 1 , N 
FF.Q=FREQ (L ) 

IF (FRQ.LE . 25. ) GO TO 1 
XX (L) =XX(L)/28. 

GO TO 8 

1 IF (FRQ. LS . 15. ) GO TO 2 

XX(L) =XX (L)/(1 05.5-3. 14*FRQ) 

YY ( L) =YY (L)/(181.32-7.588*FRQ) 

Z (L)=Z (L) / (17 7.26-7. 4 84*FRQ) 

GO TO 8 

2 IF (FRQ. LS. 10. ) GO TO 3 
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XX (L)=XX (L)/(5. 958 *FR Q-30 . 97) 

YY (L) =YY (L)/(7. 166*FRQ-39. 99} 

Z (L) =Z (L) / (6 . 4 9*FRQ- 3 2. 35) 

GO TO 8 

3 IF (FRQ. LE . 7. 5) GO TO 4 

XX (L) =XX (L)/(3.492*FRQ-6. 3 1) 

YY (L) =YY(L)/(4.252*FRQ-10. 35) 

Z (L) =Z (L) / (4. 044*FRQ- 7.89) 

GO TO 3 

4 IF (FRQ. LZ. 5.) GO TO 5 

XX (L)=XX (L)/(2.631 1*FRQ+0. 1466 7) 

YY (L) =YY (L)/(3.012*FRQ-1. 55) 

Z (L) =Z (L) / (3. 1 84*FRQ- 1.44) 

GO TO 8 

5 IF (FRQ. LE. 3.) GO TO 5 

XX (L) =XX (L)/(2 . 631 1 * F RQ + 0 . 14657 ) 

7 YY (L) =YY (L )/ ( 2. 702*FR Q) 

Z (I) =Z (L) / (2. 9 2*FRQ) 

GO TO 8 

6 XX (L) =XX(L)/(2.72*FRQ) 

GO TO 7 

8 CONTINUE 

9 CONTINUE 

C THIS BLOCK, ENDING WITH 30, THE MAGNITUDE OF THE 

C FREQUENCY DOMAIN DATA IS TAKEN AND PLACED IN A REAL 

C ARRAY, ’ZX’ FOR THE X-COIL^ZY' FOR THE Y-COIL AND ZZ 

C FOR THE Z -COIL. THIS DATA IS DIVIDED EY THE NUMBER OF 

C SAMPLE POINTS AND SQUARED TO DETERMINE POWER. 

DO 30 11= 1 ,N 

ZX (II ) =ZX (II) + (CA3S (XX (II) ) /FN| **2 

ZY (II) =ZY (II) + (CABS (Y Y (II ) ) /FN| **2 

ZZ (II) =ZZ (II) + (CAB S ( Z (II) )/FN) **2 

SX7 (I 1) =SX Y(I 1 ) + ((XX (II) /FN) *C0NJG (YY (II) /FN) ) 

COXY (II) =COXY (II) +CAB S ( (XX (II ) /FN) *CONJG (YY (II ) /FN) ) 
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SYZ (I 1) =SY Z(1 1 ) + ( (YY ( II) /FN) *CONJG (Z (II) /FN) ) 

COYZ (II) =CCYZ (II) + CA3 S ( (Y Y (II ) / FN) *CONJG (Z (I1)/FN)) 

SZX (II) =SZX(I 1) +( (Z (I 1) /FN) * C 0 N J G (XX (II) /FN) ) 

COZX (II) =COZX (II) +CA8 S ( (Z (II) /FN) *CONJG (XX (II) /FN) ) 

30 CONTINUE 
70 CONTINUE 

C THIS BLOCK AVERAGES THE DATA POINTS ADDED IN BLOCK 30 

C ABOVE BY 'NR'. AT THIS POINT THE POWER CURVES ARE 

C CONVERTED INTO POWER SPECTRAL DENSITIES. 

C THE COHERENCES AND STOKES PARAMETERS ARE 

C COMPUTED IN ALL THREE PLANES. 

DO 33 13= 1 ,N 

SEXY (13) = (REAL (SXY (13 ) ) /FNR) *4. 

S3XY (13) = (AIM AG (SXY (I 3) ) /FNR) *4 . 

S2YZ (13) = ( REAL (SYZ (13 ) ) /FNR) *4. 

S 3 YZ (13) = (AIM AG (SYZ (I 3) ) /FNR) *4 . 

S2ZX (13) = (PEAL (SZX (13 ) ) /FNR) * 4 . 

S3ZX (13) = ( AIM AG (SZX (I 3) ) / FNR) *4 . 

COXY (13 ) = (COXY (13) ) /3QRT(ZX (13) *ZY (13) ) 

COYZ ,( 13) = (COYZ (13) ) /SQRT( ZZ (1 3) *ZY (13) ) 

COZX (13) = (COZX (13) ) /3QRT(ZX(I3) *ZZ (13) ) 

ZX (13) =ZX (13) *T/FNR 
ZY (13) =ZY ( 13) *T/FNR 
ZZ (13) =ZZ ( 13) *T/FNR 
S1XY (13) = ( (ZX (13) - ZY (13)) /T) *2. 

S 1YZ (13) = ( (ZY (13) -ZZ ( 13) ) /T) *2. 

S1ZX (13 )= ( (ZZ (13) - ZX (13) ) /T) *2. 

SOXY (13) = ( (ZX (13) 4 ZY (13)) /T) *2. 

SOYZ (13 )= ( (ZZ (13) +Z7 (13) | /T) * 2. 

SOZX (13) = ( (ZX (13) fZZ(l3)) /T) *2. 

33 CONTINUE 

C THE NEXT DO LOOP CONVERTS THE B- FIELD POWER SPECTRAL 

C DENSITY INTO DECIBELS REFERENCED TO 1 NANOTESLA **2 

C PER HERTZ . 
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DO 32 1 = 1 , N 

ZX (I) =10. * ALOG 10(ZX(I ) ) 

ZY (I) =10. * ALOG 1 0 (Z Y ( I ) ) 

ZZ (I) =10. * ALOG 10 (ZZ(I ) ) 

32 CONTINUE 

C THE NEXT 3 LOCK ELIMINATES THE FIRST 10 DATA POINTS SO 

C THAT THE FREQUENCY IS .04 TO 13 HZ. 

C S1,S2, AND S3 ARE NORMALIZED BY SO. 

C DEGREE OF POLARIZATION, ORIENTATION ANGLE ALPHA, 

C AND THE E L LIPTICIT Y ARE ALSO COMPUTED. 

DC 950 14= 1,3 1 32 
FREQ (14) = F REQ (14+10) 

FRQ2 (14) =F PQ2 (14+10) 

ZX (14) =ZX (14+10) 

ZY (14) =ZY (14+ 1 0) 

ZZ (14) =ZZ (14+1 0) 

COXY (14) = C CX Y (14+ 1 0) 

SHY (I4)=S 1XY (14+10) / SOXY (14+10) 

S 2XY ( 14 ) =S2XY (14+10) / SOXY (14 + 10) 

5 3XY.(I4)=S3XY (14+10) / SOXY (14+10) 

PXY (14) =SQ RT (S 1 XY (14 ) **2+S2XY (14) **2 + S3XY (14 ) **2) 

DEXY (14) = 0 .5*ARSIN (5 3 XT (14) /PXY (14 ) ) 

SOXY (14) = 1 0. *A LOG 10(5 OXY (14+101 ) 

CO YZ (14) =C OYZ (14+10) 

S1YZ (14) =S 1YZ (14+10) / SO YZ (14+10) 

S 2YZ (14) =S2YZ (14+1 0) / SO YZ ( 14 + 1 0 ) 

S3YZ (14) = S 3YZ (14+10) /SOYZ (14+10) 

?YZ (14) =SQ RT ( S 1 YZ ( 14 ) **2+S2YZ (14) **2+S3YZ(I4)**2) 

DEYZ ( 14 ) = 0 .5* A RSIN (3 3 YZ (1 4 ) /? YZ (14) ) 

SOYZ (14) = 1 0. * A LOG 10 (SOYZ (I 4+ 10) ) 

COZX (14) =COZX (14+1 0) 

S1ZX (14 ) = S 1ZX (14+10) / SOZX(I4+10) 

S2ZX (14) =S 2ZX (14+10) / SOZX (14+10) 

S3ZX (14) = S 3ZX (14+10) / SOZX (14+10) 
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PZX (14) =S Q ET { S 1ZX (14) ** 2+ S2ZX (14) **2+S3ZX (14) **2) 
DEZX (14) = 0 .5* A ESIN (3 3 ZX (14) /PZX (14) ) 

SOZX (14) =1 0. *ALOSl0(S0ZX(I4+10) ) 

AX I (14) =AT AN (S2XY (14) /S1XY (14) t /2 . 

AYZ (14) =A T AN (S 2YZ (14 ) /S1YZ (14) ) /2. 

AZX (14) =AT AN (S2ZX(I4) /S 1Z X (14) | /2 . 

950 CONTINUE 
C 

C YES S A TEC PLOTS 
C 

N?TS= 10./D ELTAF +1 . 

C 1 NPTS ' DETERMINES NUMBER OF POINTS NECESSARY 

C FOR THE 0 TO 10 HERTZ RANGE TO BE PLOTTED. 

C FOR THE ’ITS 1 AND 1 RTB' VALUES REVIEW THE WRITE-UP 

C FOR THE SUBROUTINE PROCEDURE 'DRAW? 1 . 

ITB (3) =7 
ITB (4) =5 
ITB (7) =1 
ITB (12) =1 
RTB (1) =0. 0 
RTB (2 ) =0. 0 
RTB (3) =AL A 3(1) 

READ (5, 3000) TITLE 

CALL DRAW? (NPTS, FRQ2, ZX r I T3 r RT3) 

RTE (3) =AL A 3 (2) 

READ (5, 30 00) TITLE 

CALL DRAW? (NPTS , FRQ2, ZY ,1 TB, RTB) 

RTB (3) =AL A 3(3) 

READ (5, 3000) TITLE 

CALL DRAW? (NPTS, FRQ2, ZZ,I TB,RT3) 

RTB (3) =ALA E(4) 

READ (5, 30 0 0) TI TLE 

CALL DRAW? (NPTS, FRQ2, COXY ,IT3, RTB) 

RTB (3) =ALA 3(5) 
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READ (5, 3000) TITLE 

CALL DRAWP (N?TS,FRQ2, SOXY , ITB, RIB) 
RTB (3) =ALA B(6) 

READ (5, 3000) TITLE 

CALL DRAW? (N?TS,FRQ2, Si XY , ITB , R T 3) 
RTB (3 ) = AL A B (7) 

READ (5, 30 0 0) TITLE 

CALL DRAW? (N?TS, FRQ2, S2XY , ITB, RTB) 
RTB (3) =ALA 3(3) 

READ (5, 3000) TITLE 

CALL DRAW? (NPTS ,FRQ2 , S3 XY ,IT3 , RTE) 
RTB (3) = ALA 5(4) 

READ (5 r 30 0 C) TITLE 

CALL DRAW? (N?T S, FRQ2, COYZ ,IT3, RTB) 
RTB (3) =ALAB(5) 

READ (5, 3000) TITLE 

CALL DRAW? (N?TS,FRQ2, SO YZ, ITB, RIB) 
RTB (3) =AL A B(6) 

READ (5, 30 0 0) TITLE 

CALL DRAW? (N?TS,FRQ2, S 1 YZ , ITB, RIB) 
RTE (3) =ALA 3(7) 

READ (5, 3000) TITLE 

CALL DRAW? (NPTS ,FRQ2 , S2 YZ , ITB , R 1 3 ) 
RTE (3 ) = AL A B (8) 

READ (5, 3000) TITLE 

CALL DRAW? (N?TS , FRQ2 , S3YZ , ITB, RTB) 
RTE (3) =AL A 3 (4 ) 

READ (5, 3000) TITLE 

CALL DRAW? (NPTS ,FRQ2 , COZX , ITB, RTB) 
RTB (3 ) = AL A E (5) 

READ (5, 30 0 0) TITLE 

CALL DRAW? (NPTS, FRQ2, SOZX , ITB, RTB) 
RTB (3) =AL A B (6) 

READ (5,3000) TITLE 
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CALL DRAW P (NPTS, FRQ2, S1ZX, ITB, R 13) 
RTE (3 ) = AL A B (7) 

READ (5, 30 0 0) TI TLE 

CALL DRAW? (NPTS,FRQ2, S2ZX , ITB, RTS) 
RTB (3) =ALA 3(8) 

READ (5, 3000) TITLE 

CALL DRAW? (N?TS,FRQ2, S3ZX ,IT3,RT3) 
P.T 5 (3 ) = AL A E (9 ) 

READ (5, 30 0 0) TITLE 

CALL DRAWP (NPT S , FRQ2 , AX Y, ITB, R 1 5) 
RTE (3) =AL A B ( 1 0 ) 

READ (5, 3000) TITLE 

CALL DRAW? (NPTS ,FRQ2 , ?XY,IT3,RT3) 
RTE (3 ) = AL A E (1 1) 

RE AD (5, 30 0 0) TITLE 

CALL DRAWP (NPTS, FF.Q2, DEXY , ITB , S TB) 
RTE (3) =ALA 5(9) 

READ (5, 3000) TITLE 

CALL DRAW? (NPT S ,FRQ2 , A YZ r I TE , 3 T 3) 
RTS (3) =ALAE(10) 

READ (5, 30 0 0) TITLE 

CALL DRAW? (NPT S , FRQ2 , PYZ,IT3,RT3) 
PTE (3) = AL A B (1 1) 

READ(5,300 0) TITLE 

CALL DRAWP (NPT S , FRQ2 , DE YZ ,ITB, RTB) 
RTB (3) = ALA 3(9) 

READ (5, 3000) TITLE 

CALL DRAWP (NPTS , FRQ2 , AZX, ITB, RTB) 
RTB (3) =ALA B (1 0) 

READ (5,3000) TITLE 

CALL DRAW? (N?TS,FRQ2 , PZX,IT3,RT3) 
RTB (3) = ALA E (11) 

READ (5, 30 00) TITLE 

CALL DRAW? (NPTS, FRQ2 , DEZX ,IT3,3TE) 



79 



3300 FORMAT (6A8 ) 
STOP 
END 



//LKED. SYSUT1 DD S?ACE=(1024, (200,203) ) 



//SO. SYS IN DD * 




LA MESA SITE, 16 SEP 82, 


1250-1330 LOCAL 


PSD OF X-COIL. AMP IN D5 


REF NT**2/HZ VS HZ 


LA MESA SITE, 16 SEP 82, 


1250-1330 LOCAL 


PSD OF Y-COIL. AMP IN D8 


REF NT**2/HZ VS HZ 


LA MESA SITS, 16 SEP 82, 


1250-1330 LOCAL 


PSD OF Z-COIL. AMP IN DB 


REF NT* *2 /HZ VS HZ 


LA MESA SITE, 16 SEP 32, 


1250-1330 LOCAL 


COHERENCE OF X S Y COILS 


VS FREQ (HZ) 


LA MESA SITE, 16 SEP 82, 


1250-1330 LOCAL 


S0( 1 0 LOG) VS FRSQ(HZ) 


X-Y PLANS 


LA MESA SITE, 16 SEP 82, 


1250-1330 LOCAL 


SI/SO VS FREQ (HZ) 


X-Y PLANE 


LA MESA SITE, 16 SEP 32, 


1250-1330 LOCAL 


32/S0 VS FREQ (HZ) 


X-Y PLANS 


LA MESA SITE, 16 SEP 82, 


1250-1330 LOCAL 


S3/S0 VS FREQ (HZ) 


X-Y PLANE 


LA MESA SITE, 16 SEP 82, 


1250-1330 LOCAL 


COHERENCE OF Y & Z COILS 


VS FREQ (HZ) 


LA MESA SITE, 16 SEP 82, 


1250-1330 LOCAL 


SO ( 1 0 LOG) VS FREQ (HZ) 


Y-Z PLANE 


LA MESA SITS, 16 SEP 82, 


1250-1330 LOCAL 


SI/SO VS FRSQ(HZ) 


Y-Z PLANE 


LA MESA SITE, 16 SEP 82, 


1250-1330 LOCAL 


S2/S0 VS FREQ (HZ) 


Y-Z PLANE 


LA MESA SITE, 16 SEP 82, 


1250-1330 LOCAL 


S3/S0 VS FREQ (HZ) 


Y-Z PLANE 


LA MESA SITS, 16 SEP 82, 


1250-1330 LOCAL 


COHERENCE OF Z 5 X COILS 


VS FREQ (HZ) 


LA MESA SITS, 16 SEP 82, 


1250-1330 LOCAL 
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SO ( 1 0 LOG ) VS FRSQ(HZ) Z-X PLANE 

LA MESA SITE, 16 SEP 82, 1250-1330 LOCAL 
SI/SO VS FREQ (HZ) Z-X PLANE 

LA MESA SITE, 16 SEP 82, 1250-1330 LOCAL 

S2/S0 VS FREQ (HZ) Z-X PLANE 

LA MESA SITE, 16 SS? 82, 1250-1330 LOCAL 
S3/S0 VS FREQ (HZ) Z-X PLANS 

LA MESA SITE, 16 SEP 82, 1250-1330 LOCAL 
A XI VS FREQ (HZ) 

LA MESA SITE, 16 SEP 82, 1250-1330 LOCAL 

DEGREE OF POLARIZATION VS FREQ(HS) X-Y 
LA MESA SITE, 16 SEP 82, 1250-1330 LOCAL 
ELL ITICITY VS FREQ (HZ) X-Y PLANS 

LA MESA SITE, 16 SEP 32, 1250-1330 LOCAL 

AYZ VS FREQ (HZ) 

LA MESA SITE, 16 SEP 82, 1250-1330 LOCAL 
DEGREE OF POLARIZATION VS F ESQ (H Z) Y-Z 
LA MESA SITE, 16 SEP 82, 1250-1330 LOCAL 
ELL ITICITY VS FREQ (HZ) Y-Z PLANS 

LA MESA SITE, 16 SEP 82, 1250-1330 LOCAL 
ASX VS FREQ (HZ) 

LA MESA SITE, 16 SEP 82, 1250-1330 LOCAL 

DEGREE OF POLARIZATION VS FREQ(HZ) Z-X 
LA MESA SITE, 16 SEP 82, 1250-1330 LOCAL 

ELL ITICITY VS FREQ (HZ) Z-X PLANE 

/* 

//GO. FT20F001 DD UN IT=3 4 00 - 4 , VO L=SER = 3 MDT5 2 , DIS P= (OL D , KEEP) , 
// LAESL= (1 ,NL, , IN) , 

// DC E= ( RSCFM=F3 ,LRSCL = 32, BLKSIZE=51 2 , DEN = 2) 

/* 

// 
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APPENDIX B 



CHEW^S RI DOE Sir E ?R2?ARATI0N 

Chew's Ridge is the site of previous land-side 
geomagnetic data collection. It is located approximately 
twenty-three air miles southeast of the Naval Postgraduate 
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school. A detailed description of the physical layout and 
operation of the site can be found in [Ref. 6]. 



The original intention of this research was to collect 
data from this remote site to provide the simultaneous data 
for the sea trials. Preparation of the site and the 
necessary equipment continued through most of the project. 
The da~a transmission from this site to the Naval 
Postgraduate School receiving location proved impossible 
with the transmitter/receivers available. 

Previous data collected from Chew's Ridge was collected 
using a Motorola system that suffered from overheating 
failures when required to operate in a continuous mode. 
Ir.tiallv another attempt was made to repair the Motorola 
unit but after one last failure this proved economically 
impossible and it was felt that this unit did not have the 
bandwidth to accomadate the addition of the PCM to the 
system. Consequently a RSPCO Model BIO transmitter-receiver 
pair was ordered. These on arrived the first of September 
and initial testing showed that, although they were capable 
of operating for long periods of time, they did not have the 
required 20 KHz bandwidth to accomadate the PCM data. 
Efforts are currently underway to expand the bandwidth to 
meet the requirements. 
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To properly orient the sensing coils at this site it is 
necessary to have an accurate known direction. In order to 
accomplish this, inquiries were mads to the National Forest 
Service to determine if tksy have survey data for the area. 
Whan it proved nonexistent, a survey by resection was 
conducted. A theodolite was used to measure relative angles 
to several prominent terrain features in the area from the 
Chew’s Ridge Site. These measurements were made three times 
and then averaged. The average relative angles were then 
drawn or. a piece of overlay paper and oriented on a 1:50,000 
map of the area. This allowed for tie best fit of the data 
and prevented blunders in the recognition of the terrain 
features. The most prominent terrain feature (Uncle Sam 
Mountain) was then chosen for the end of the orienting line. 
The Grid Azimuth to Uncle Sam Mtn. was then computed from 
the coordinates of the two locations. Finally the True 
Azimuth was calculated to be 257.4 degrees. Estimated 
• accuracy of this measurement is one degree. 
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